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, Abstract 

£f) ■ We present a derivation of Wandzura-Wilczek (WW) like relations for skewed par- 

ton distributions. It is demonstrated for photon-pion scattering that the skewed twist-3 
parton distributions contributing to the DVCS amplitude have discontinuities at the 
points x = ±£ in the WW approximation. This may lead to a violation of factorisation 
for the twist-3 DVCS amplitude with transverse polarization of the virtual photon. 
Qh! We show, however, that the contribution of the divergencies to the scattering of a 

Q_ ( ' transversely polarized virtual photon affects DVCS observables only at order 1/Q 2 and 

can be neglected at twist-3 accuracy. For the scattering of a longitudinally polarized 
photon the twist-3 amplitude is free of such divergencies. 
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Introduction 

Deeply Virtual Compton Scattering (DVCS) [[|, || is the cleanest hard process which is 
sensitive to the Skewed Parton Distributions (SPD) and has been the subject of extensive 
theoretical investigations for a few years. First experimental data became recently available 
(see e.g. J3], |J) and much more data are expected from JLAB, DESY, and CERN in the 
near future. It was demonstrated [| [jj that in leading order (1/Q°) the DVCS amplitude 
factorizes. However, as the typical Q 2 are by no means large, studies of the power corrections 
to the DVCS amplitude are very important. 



Recently the DVCS amplitude was computed [§, |9|, |10[ including the terms of order 
0(1/Q). The inclusion of such terms is mandatory to conserve the electromagnetic gauge 
invariance of the DVCS amplitude at this order of 1/Q. Moreover, they provide the leading 
contribution to some spin asymmetries. To this order the amplitude depends on a set of 
new skewed parton distributions. Recently it was shown in ref. |10| that in the so-called 



Wandzura-Wilczek (WW) approximation these new functions can be expressed in terms of 
twist-2 SPD's. WW relations for SPD's were also discussed in ref. |]TTJj , however the authors 
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of this paper did not take into account operators which are total derivatives. Obviously 
these operators are crucial for the description of non-forward matrix elements. 

In this note we give a derivation of the WW-like relations which is technically slightly 
different from the approach used in ||10|| , see the Appendix for details. Then we demonstrate 
for the case of DVCS on the pion the new SPD's in WW approximation generically posses 
discontinuities at the points x = ±£. This leads to a formally divergent expression for the 
DVCS amplitude, which however contributes to the DVCS observables at the accuracy l/Q 2 
and hence beyond our accuracy. 

We also observed that in the twist-3 DVCS amplitude^ with longitudinally polarized 
virtual photons the divergencies mentioned above are cancelled. We discuss the theoretical 
implications of this phenomenon and its possible experimental verification. 



DVCS on the pion 



The expression for the DVCS amplitude on the pion as obtained in || and reproduced in 
refs. H [UJ has the form: 
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where P = (p + p')/2, Q = (q + q')/2, A = q — q' = p' — p = —2t;P + A±, while p, p', q, q' are 
the initial and final momenta of pion and photon, respectively. The three terms in the square 
brackets are individually gauge invariant up to order A^_, i.e. q^T^ = 0(A 2 ± ) = q^T^ v . The 
second and third term contain the new twist-3 contributions to the DVCS amplitude. As P 
and Q are longitudinal and transverse the second corresponds to longitudinally polarized 
virtual photon and third to transverse polarization. 
The amplitude depends on new SPD's defined asf] 
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Here the light-cone vector n is normalized as n ■ P = 1 



Recently in ref. [TI| it was suggested to use the Wandzura-Wilczek approximation to 
express functions like H 3 in terms of twist-2 function Hi. This approximation gives the 
kinematical part of twist 3, while the genuine, dynamical, higher twist is neglected. If 



*We use the notion of 'twist-3 amplitude' as synonymous to 'amplitude of order l/Q' 
t Note that the function Ha in ref. is defined with opposite sign 
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one adopts the operator definition of twist, this approximation, because of neglecting the 
operators containing the gluon field strength, corresponds to twist 2. In the Appendix we 
give alternative derivation of such relations and in the next section we show that the WW 
approximation generically leads to formally divergent expression for the DVCS amplitude. 
We also specify the twist-3 helicity amplitudes which are free of such divergencies. 

Discontinuities of the twist-3 distributions in WW approximation 

In this section we study the properties of the twist-3 distributions in the WW approximation. 
We shall see that generically they posses discontinuities at the points x = ±£. We shall 
illustrate this by the example of H 3 (x,£). 

Using the general expression ( |A.11|) discussed in the Appendix it is easy to obtain the 
function Hs(x,£) in WW approximation. It can be expressed in terms of the twist- 2 SPD 
Hi(x,£) as follows (c.f. |T0| ): 
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Note that in the limit £ = this reduces to an expression, which is similar to the standard 
WW approximation for the transverse spin structure function gx- Note that the specific form 
of (4) guarantees the correct symmetry properties of H3 || , which is a sensitive check of our 
result. Another interesting limit to discuss is £ = 1. In this limit our expression is similar 
to the corresponding relation for the light cone amplitude for vector mesons (which is the 
closest case worked out in literature, namely by Ball and Braun [H]]). The relation (H) also 
allow to generalize the WW relation for vector mesons distribution amplitudes obtained in 



L7\[ to the case of a meson of arbitrary spin. This can be done with help of crossing relations 
between pion SPD's and distribution amplitudes of resonances |14|l - 

From eq. (|^) one can derive the following relations for the Mellin moments of Hi and 

# 3 : 



j\xH 3 (x,0 = -^J\xH 1 (x,O = 



(5) 



where the last equation follows from the polynomiality condition for SPDs . This relation 
is an analog of the Burkhardt-Cottingham sum rule, it is valid beyond WW approximation. 
This particular moment, like the other even moments, corresponds to the non-singlet part 
of the SPDs, while the contribution of the singlet part is zero for the trivial reason that 
the pion SPD antisymmetric in x. Another interesting relation, which is an analog of the 
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Efremov-Leader-Teryaev sum rule [ IB] , can be obtained for the second Mellin moments of 
the SPD's, corresponding to the singlet part, appearing in DVCS: 

rl i d rl i 

J ^dxx H 3 (x,£) = -- — J^dxx #i(x,0 = -£M 2 , (6) 

where M 2 is a momentum fraction carried by the quark in the pion. To obtain the last 
equality we used the generalized momentum sum rule for SPDs at A 2 = as derived in ITS 



Let us stress that this relation is valid beyond the WW approximation. Its generalizations 
to the nucleon case was recently discussed in 

With the help of expression (|4]) we can easily compute the behaviour of H3 near the 
points x = ±£. Let us consider the difference of left and right limits of the function H 3 at 
x — > ±£, the result is: 
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where vp J denotes a principle value integral (valeur principal). From these expressions we 
see that for a wide class of functional forms for H\ the corresponding function H3 exhibits 
discontinuities at the points x = ±£. Even for smooth functions Hi the discontinuities are 
non-zero. To clarify this point, it is instructive to integrate by part, imposing the natural 
condition Hx(l, £) = #i(-l, f) = 0. One gets: 

lim{<^ + *,fl -<«-(<:-«)} = J-± vp f J uH ^\ (9 ) 

hm{<-(-« + «)-«,-(-«-«)} = (10) 

The sum of the two jumps which arises in the convolution integral for the amplitude can be 
expressed in terms of the real part of the twist- 2 DVCS amplitude as follows: 

hm{# 3 W (£ + 5, - Hf w {i - 5, + Hf w {-$ + 5, - Hf w {-t -6,$} 

= -~ vp f duH^oi^—, + -^-7}. (11) 

The latter expression is nothing else than the derivative of the real part of twist two DVCS 
amplitude. As DVCS is an observable physical process (and as already been observed) the 



integral on the right hand side of (11) cannot be zero and because the physics has to show 
a non-trivial dependence on the skewedness parameter it cannot be independent of £. This 
shows that the r.h.s. of ( pTTj ) cannot vanish and therefore H 3 really has to show discontinuities. 
The expression for the WW- part Hj[ £) of the function Ha(x,£) is slightly different: 
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It can be checked that the function H A vw (x, £) also posses discontinuities at the points 
x = ±£. The expression for these discontinuities has the form: 



hm{i^(±£ + 5,0- Hj w {±t - 5, 



1 d f 1 
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Note that the difference of these discontinuities (which is the combination entering the 
amplitude (|I|)) coincides with (|TTj). The presence of discontinuities in the functions H 3 and 
Ha exactly at the points x = ±£ leads to a formally divergent result for the amplitude at 
the order 0(1/Q), see expression eq. (|l|). This indicates a violation of factorisation in order 
0(1/ Q) for the DVCS amplitude in WW approximation. The considered violation disappears 
for virtual final photon, as the pole of the quark propagator no more occurs at the point 
x — £. Thus the production of Drell-Yan pairs does not pose any problems. Experimental 
studies comparing real and virtual photon production thus provide an excellent opportunity 
to check factorisation. 

It is important to realize that the formal divergencies of the amplitude in the WW 
approximation are canceled for certain combinations of helicity amplitudes because the jumps 
in H 3 and Ha are related to each other, see eqs. (0,0). Since the divergencies occur 
only at the point x = ±£ it is sufficient that the jumps cancel at this specific value to 
save factorisability. Specifically there is no problem for the amplitudes with longitudinal 
polarization of the virtual photon, because the corresponding amplitude is proportional to 
the following combination of the functions H 3 and Ha'- 
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Thus the discontinuities generate divergencies only for the twist-3 DVCS amplitude with 
transversely polarized virtual photon. The contribution of the "problematic" part of the 
twist-3 DVCS amplitude (corresponding to scattering of the transversely polarized virtual 
photons) to the differential cross section is suppressed after the contraction with the polariza- 
tion vector of the emitted real photon by two powers of the hard scale, i.e. by l/Q 2 , relative 
to the leading order result and do not contribute to observables in order 0(1/ Q). The DVCS 
differential cross section to this accuracy gets contributions only from the longitudinal part 
of the twist-3 amplitude which is free from divergencies. 

The physical reason of these peculiarities for the transverse photon case is rather similar to 
the origin of the famous Callan-Gross relation. The dangerous pole of the quark propagator, 
which in combination with jumps lead to the violation of factorisation, corresponds to an 
on-shell quark, which may absorb only transverse photons. We therefore expect that this 
situation should persist also in the case of nucleon and deuteron targets. This would be 
interesting to check by an explicit calculations. 
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The considered relationship between twist-2 and twist-3 terms can also be discussed with 
regard to electromagnetic gauge invariance. The leading twist-2 term is gauge invariant 
only if one neglects the transverse component of the momentum transfer Aj_. In 0(A±) 
accuracy it requires the consideration of a quark-gluon diagram, whose contribution by use 
of the equation of motions is expressed through the twist-2 SPD Hi and the new twist three 
SPDs H 3 , Ha. While the Hi term combines with the leading result to a gauge- invariant 
expression (see two first lines in eq. (p])), as anticipated in ref. ||13|| , the other terms provide the 
additional contribution which violates factorisation in the WW approximation for transverse 
polarization of the virtual photon. 

Generally, the appearance of the jumps just at the points x = |£| is not unnatural from 
the physical point of view, as this is just a transition point between the regions, where SPD 
has quite different physical meaning || , accommodated, in particular, in the two-component 
model of SPDs JRJ. 

As final remark we note that preliminary estimates of the function H 3 at a low normaliza- 
tion point in the instanton model give a function without discontinuities and hence indicates 
that the WW approximation for SPD's is not valid within this specific model. 

Conclusions 

We demonstrated by explicit calculations that Wandzura-Wilczek like relations for SPD's 
entering the description of the DVCS amplitude at order 0(1/ Q) lead to a violation of 
factorisation for the twist-3 DVCS amplitude with transversely polarized virtual photons. 
However, one can easily see that the dangerous divergencies do not contribute to DVCS 
observables at the order 1/Q but at the order 1/Q 2 . Therefore these divergencies affect only 
twist-4 corrections which are beyond the scope of the present paper. One cannot exclude 
that the kinematical contributions of twist 4 will cancel the considered divergencies. This 
promising opportunity is suggested by the paper p0| |, in which a part of the 1/Q 2 term is 



identified which makes the contribution of the jumps equal to zero after the contraction with 
the real photon polarization vector is performed. 

The divergencies in the amplitude are related to the fact that the additional functions 
obtained with the help of WW-like relations generically contain discontinuities at the points 
x = ±£ what in turn leads to formally divergent results for the DVCS amplitude with 
a transversely polarized photon. However the divergencies contribute to the observables 
only at the accuracy 1/Q 2 and can be neglected at twist-3 accuracy. We observed that the 
divergencies are canceled in the twist-3 DVCS amplitude with a longitudinally polarized 
virtual photons. Such a cancelation of the divergencies in the longitudinal twist-3 amplitude 
implies that the DVCS observables up to accuracy 0(1/Q) can be estimated using the WW 
approximation, at least at leading order in a s . The case of DVCS off the nucleon will be 
studied elsewhere. 
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APPENDIX: Derivation of the WW relations 



Here we briefly describe the method used to derive the WW relations for skewed parton 
distributions. Our approach is very close to approach used in The difference lies only 
in technical details. We have also been informed that similar results have been obtained 
independenly in [p0|. 



Our starting point is the following equations derived in [19 
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where we do not write out explicitly all the path-ordered gauge factors, and V a denotes the 
total derivative defined as : 
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with generic Dirac matrix structure T and [x, y\ = Pexp [igj^dt (x—y) fl A^(tx + (l — t)y)]. Note 
that in the matrix elements the total derivative can be easily converted to the momentum 
transfer: 



(p'\V^ifj(tx)r[tx, —tx]i[)(—tx)\p) = i(p' — p)^{p'\ip{tx)T[tx, —tx]i/j(—tx)\p) 



(A.4) 



The general method to obtain Wandzura-Wilczek relations is to take the matrix elements 
in the LHS and RHS ( |A.l ) and ( |A.2| ) and insert their parametrisation. This provides us 
with a system of equations for the twist three functions like H 3 and Ha- Such a method 
was effectively also used for the p-meson distribution amplitudes, see |0| [L8| . In the case of 



skewed distributions it is more convenient to solve the system of equations ( |A.1| ), ( |A.2 ) at 
the operator level and then take the matrix elements. Operator solution means that one has 
to express non-symmmetrical operator %p{x)^( ^ b ip{—x) through the two point symmetrical 
operators %[)(x)tf{^z)ip(—x) and three point quark-gluon operators. 
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Consider as an example the solution for the vector operator. Substituting ( |A.2| ) into 
( |A.1| ) we obtain an equation which contains only one non-symmetrical vector operator : 



ui){ux)-i^{-ux) = [{xVf - x 2 {V) 2 } j dtt(u-t)i[j(tx)^ fl ip(-tx) + 
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dx k Jo 

where ellipses stands for the contributions of the three point quark-gluon operators. It is 
convenient to rewrite this equation in the compact form: 

f(u) = k 2 / (u-t)f(t) + <p(u) (A.6) 



where we introduced the following notations: 

f(u) = ui>(ux)-f^i/j(-ux), /(0) = 0, 



k = ^J(xV) 2 -x 2 (V) 2 (A.7) 

and (p(u) denotes the known terms in the RHS of the (|A.5|) . Our aim is to solve this integral 
equation and find f{u). It it easy to see that ( |A.6| ) can be reduced to the simple second 
order differential equation: 

f"(u) = k 2 f(u) + V "(u), (A.8) 

with the boundary conditions: 

/(0) = 0, f(0) = A^(o)^(o) (A.9) 



The solution is 
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Substituting the definitions (|A.7|) we obtain an explicit expression. We then simplify the full 
answer by neglecting the square of the momentum transfer. This corresponds to the formal 
substitutions T> 2 = 0, k — xT>. Then after simple manipulations we obtain the following 
operator equation: 



\ f da \aV, (e-^ xV) - e^ xV) ) + (e~^ xV) + e* xV) ) 
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(xV) _ a(xV)\ 
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Jo Jo ox k 

where at = 1 — a and ellipses are contributions of the three point quark gluon operators 
which we do not write explicitly for the sake of simplicity. This operator relation can be 
used for the derivation of various WW relations, because in the RHS of the eq. (|A.11|) only 
two-point quark operators of twist- 2 contribute in this case. Also the general solution ( |A.10[ ) 
can be applied for all-order summation of the kinematical power suppressed terms in the 
DVCS amplitude. 
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